Renormalization group study of random quantum magnets 
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We have developed a very efficient numerical algorithm of the strong disorder renormalization 
group method to study the critical behaviour of the random transverse-field Ising model, which is a 
prototype of random quantum magnets. With this algorithm we can renormalize an TV-site cluster 
within a time iVlogiV, independently of the topology of the graph and we went up to N ~ 4 x 10 s . 
We have studied regular lattices with dimension D < 4 as well as Erdos-Renyi random graphs, 
which are infinite dimensional objects. In all cases the quantum critical behaviour is found to be 
controlled by an infinite disorder fixed point, in which disorder plays a dominant role over quantum 
fluctuations. As a consequence the renormalization procedure as well as the obtained critical prop- 
erties are asymptotically exact for large systems. We have also studied Griffiths singularities in the 
paramagnetic and the ferromagnetic phases and generalized the numerical algorithm for another 
random quantum systems. 
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I. INTRODUCTION 

Quantum phase transitions take place at T = tem- 
perature by varying a control parameter, which is in- 
volved in the Hamiltonian of the system^. Experi- 
mental examples in which quantum phase transitions 
play an important role are among others rare-earth 
magnetic insulators^, heavy- fermion compounds^, high- 
temperature superconductors^ and two-dimensional 
electron gases^. Generally, when quantum fluctuations 
are weak the ground-state of these systems is ordered, 
whereas for strong quantum fluctuations we are in the 
quantum disordered region. In between there is a quan- 
tum phase-transition point, the effect of which is mani- 
fested also in finite temperature in the so called quantum 
critical region. 

A paradigmatic model having a quantum phase- 
transition is the Ising model in the presence of a trans- 
verse field of strength, h, where h plays the role of the 
control parameter. For h < h c (h > h c ) the system is fer- 
romagnetic (paramagnetic) and at h = h c there is a quan- 
tum critical point. Experimentally this system is real- 
ized by the compound LiHoF 4 , which is a dipole coupled 
Ising ferromagnet, and which is placed in a magnetic field 
transverse to the Ising axis of strength, H t . This mag- 
netic field splits the ground-state doublet and therefore it 
acts as an effective transverse field of strength h ~ H? . In 
the above compound one can introduce randomness, by 
substituting the magnetic Ho by a nonmagnetic Y. Then 
the obtained system LiHo x Yi_ :c F4 is the experimental 
realization of a random quantum magnet^. However, 
the transverse field induces also a random longitudinal 
fields via the off-diagonal terms of the dipolar interac- 
tion, therefore there are several open questions both the- 
oretically and experimentally about the low-temperature 
behavior of this compound. 

Here we consider a theoretically simpler problem, the 
random transverse-field Ising model (RTIM) , which is de- 



fined by the Hamiltonian: 



H = - ^ JijVitf 



(1) 



Here the of if 



are Pauli-matrices at sites i (or j) of a 
lattice and the nearest neighbour couplings, Jy, and the 
transverse fields, hi, are independent random numbers, 
which are taken from the distributions, p(J) and q(h), 
respectively. In this paper we have used two different 
disorder distributions in order to check universality of 
the properties at the critical point. For both disorders 
the couplings are uniformly distributed: 



P (j) = e(j)e(i - j) , 



(2) 



where &(x) denotes the Heaviside step-function. 

For box-fe disorder the distribution of the transverse- 
fields is uniform: 



q(h) = —e(h)Q(h b - h) , 
h b 



(3) 



whereas for fixed- h disorder the initial values of the 
transverse-fields arc constant: 



q(h) = 6{h - h f ) 



(4) 



The quantum control parameter is defined as 9 — log(/i;,) 
and 9 = log(/i/), respectively. 

Theoretical study of the critical behaviour of the RTIM 
is a complicated issue, since one should treat the joint 
effect of quantum and disorder fluctuations, as well as 
non-trivial correlations. Most of the results in this field 
are known in one dimension (ID), due to a renormaliza- 
tion group (RG) treatmenljii, which has been introduced 
by Ma, Dasgupta and Hu^ and further developed by 
Fisher—. This RG procedure works in the energy space 
and in each step the largest local parameter in the Hamil- 
tonian, either a coupling or a transverse field is decimated 
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out and at the same time new effective parameters are 
generated perturbatively. Repeating this procedure the 
energy scale in the system, measured by the largest ef- 
fective parameter goes to zero and one studies the dis- 
tributions of the renormalized couplings and that of the 
transverse fields at the fixed point. 

In ID where the topology of the system is preserved 
during renormalization Fisher— has solved analytically 
the RG equations and in this way precise information has 
been obtained about the behavior of the RG flow in the 
vicinity of the fixed point . One important result, that the 
distribution of the renormalized parameters (couplings 
and/or transverse fields) at the fixed point are logarith- 
mically broad. This means that disorder fluctuations are 
completely dominant over quantum fluctuations and the 
perturbative treatment used in the calculation of the new 
parameters is asymptotically exact. In this, so called infi- 
nite disorder fixed pointii (IDFP) of the ID model exact 
results has been derived for the critical exponents, as well 
as for several scaling functions. 

In more complex geometries, even for a ladder with 
w > 1 legs the topology is changing during the 
RG process and therefore one relies on numerical im- 
plementations of the RG procedure. For a ladder 
with a finite width this procedure can be performed 
straightforwardly-^, during which the system is renor- 
malized to an effective chain. As a result the critical 
singularities of ladders are identical to those of a chain, 
however, from the w-depcndcncc of the amplitudes one 
can deduce cross-over functions and estimate the singu- 
larities of the 2D system through finite size scaling. 

If the sample is isotropic, say an L x L part of a 2D 
lattice one can not successfully renormalize it numeri- 
cally with a naive application of the RG rules. It is due 
to the fact that after /i-decimation steps a large num- 
ber of new effective couplings are generated and the sys- 
tem will soon look like to a fully connected cluster, for 
which the further decimations are very slow. To speed 
up the process new ideas have been introduced, one of 
those is the so called maximum rulei^. This is applied 
in a situation, when between two sites two or more cou- 
plings are generated and the maximum of those is taken. 
This procedure is certainly correct at an IDFP, where the 
couplings have typically very different magnitudes. To- 
gether with the maximum rule one can also use some fil- 
tering conditional to eliminate the latent, non-decimated 
couplings of the Hamiltonian, which are at such a posi- 
tion, where a larger parallel coupling will be generated. 
Using these tricks one could numerically study the 2D 
systemiSiiSr^i up to a linear size L = 128 — 160. The 
critical exponents calculated in this way are in agreement 
with the finite-size scaling extrapolations obtained in the 
strip geometry^ and with the results of quantum Monte 
Carlo simulations^. The correctness of the method is 
also checked by comparing results of Monte Carlo simu- 
lations about the random contact process^, which is a 
basic model of reaction-diffusion processes in a random 
environment. According to an RG analysis the critical 



behavior of this model (for strong enough disorder) is 
identical to that of the RTIM^, which is indeed found in 
ID and 2D calculations. 

The model of real random quantum magnets, the 
RTIM in 3D can not be studied even by the above re- 
fined RG algorithm, since the available finite sizes of the 
system are too small to obtain stable estimates about the 
properties of the fixed point. In early studies the possi- 
ble presence of an IDFP is expected^*, but no evidence 
in favour of this conjecture has been presented. Also no 
studies are available about the random contact process 
in 3D. In higher dimensions no results of any kind are 
known, thus it is a completely open question, if there is 
an upper critical dimension, D Ul such that for D < D u 
infinite disorder scaling works and for D > D u we have 
conventional random criticality. Furthermore, it is also 
an open question, if for D > 2 the fixed point is universal 
and does not depend on the actual form of the disorder. 

A possible way to answer to the questions presented 
in the previous paragraph is to improve the numerical 
algorithm of the RG procedure and make it capable to 
study three- and higher dimensional systems. In this pa- 
per we describe such an improved algorithm, which uses 
the maximum rule but otherwise has no further approx- 
imations. This means that the results of our algorithm 
are identical to that of any naive implementation of the 
RG method (having also the maximum rule) for any fi- 
nite graphs, say with N sites and E edges. However, we 
gain considerable time in performance: while the naive 
method works in t ~ 0(N 3 ) time, this is for the improved 
algorithm t ~ 0(N log N + E). Having this performance 
at hand we could treat finite clusters up to 4 x 10 6 sites. 

We have used this algorithm to study the properties of 
the RTIM in different dimensions. Some preliminary re- 
sults of these investigations have already been presented 
elsewhere. The critical behaviour of the 3D and 4D sys- 
tems is announced ir>2£. Here we give details of the deter- 
mination of the critical parameters as well as analyze the 
scaling behavior of different quantities in the off-critical 
region, too. The model in 2D has been studied by another 
algorithm ir>ii, which algorithm has basically the same 
performance as the present one in 2D, however, which is 
less effective in higher dimensions. Here the 2D results 
are merely used to compare those with higher dimen- 
sional results. We have also studied Erdos-Rcnyi ran- 
dom graphs^ with a finite coordination number, which 
arc infinite dimensional objects. In this way we have got 
information about the possible value of the upper critical 
dimension, D„, in our system. 

Our paper is organized as follows. The RG procedure, 
the basic decimation rules and the essence of the im- 
proved algorithm are presented in Sec|H] Results about 
the critical behaviour of the system in different dimen- 
sions are presented in Scc lIIII The results are discussed in 
Sec II VI and possible extension of the improved algorithm 
for another models are given in the Appendix. 
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II. THE RENORMALIZATION GROUP 
METHOD AND THE IMPROVED ALGORITHM 

The so called strong disorder RG (SDRG) method^ 
has been introduced by Ma and coo-workers^ to study 
ID random antifcrromagnctic Hciscnbcrg chains. For the 
RTIM it was Fisher—, who used first this method and 
solved analytically the RG equations in ID. Here we dis- 
cuss the basic steps of the RG procedure for the RTIM 
model. 

The SDRG method works in the energy space: at each 
step the largest parameter of the Hamiltonian in Eq. ([1) , 
which is denoted by D, = max{ Jy, hi}, is eliminated. 
Here we should consider two possibilities. 

./-decimation : If the largest term is a coupling, say 
fl = Jij , than the two connected sites i and j are merged 
into a spin cluster having an effective moment ft = /i; + fij 
(in the original model pn = l,Vi), which is placed in an 
effective transverse field of strength: h = hihj/Jij. This 
latter formula is obtained in second-order perturbation 
calculation. 

fe-decimation : If the largest term is a transverse field, 
say n = hi, than this site brings a negligible contribu- 
tion to the (longitudinal) susceptibility of the system and 
therefore decimated out. At the same time new cou- 
plings are generated between all sites, say j and k, which 
were nearest neighbours to i. In second-order perturba- 
tion calculation these couplings are given by JjiJik/hi. If 
there is already a coupling, Jjk > 0, between the two sites 
we use the maximum rule: Jjk = ma,x{JjiJik/hi, Jjk}- 

In the naive application of the SDRG rules for higher 
dimensional clusters there is a problem with the h- 
decimation steps, during which several new couplings are 
generated. As a result our cluster will be transformed 
soon into an almost fully connected graph, having 0(N 2 ) 
edges. Consequently at any further decimation step one 
needs to perform 0(N 2 ) operations, which leads to a 
performance in time t ~ 0(N 3 ). 

A. Improved algorithm 

Here and in the following we assume, without restrict- 
ing generality, that Jy < l,Vi,j. Then we define the 
set of local maximum, which consists of such parameters, 
which are larger (not smaller) than any of its neighbour- 
ing terms. Considering a coupling Jy is a local maxi- 
mum, provided Jy > hi, Jij > hj, and Jy > Jik, Vfc, 
as well as Jy > Jy, VZ. Similarly a transverse field, hi, 
is a local maximum, if hi > Jy, Vj. We have shownii 
that the local maximum can be decimated independently, 
the renormalization performed in any sequence gives the 
same final result. 

In the improved algorithm we concentrate on the trans- 
verse fields, the decimation of which being the most dan- 
gerous in respect of the performance of the algorithm. 
Our strategy is to avoid any Zi-decimation during the 
renormalization. For this purpose we divide the sites 



into two classes. For 'inactive' sites the transverse field 
is a local maximum and the site has a weight, Zj = 1, 
whereas all the remaining sites are termed as 'active' hav- 
ing weights, li = 0. In the log-energy space between sites 
i and j we define a distance, c!y > 0, as: 

eZy = - In J^ + ^ In h t + In hj . (5) 

Having an inactive site, k, between i and j and decimat- 
ing it out the RG rules lead to the additivity property: 
cZy = dik + dkj , which - according to the maximum rule - 
should be compared with dy in Eq.([5]) and their minimal 
value is taken. Generally, the true distance between i and 
j , denoted by #y is given by the shortest path which goes 
over the inactive sites. It is also easy to see, that deci- 
mating out all or a subset of inactive sites is equivalent 
to find in the original problem the shortest paths among 
the non-decimated sites which go through the decimated 
sites. This is a well known graph-theoretical problem^! 
for which efficients numerical algorithms are available^ 
In general, however, one should also deal with the active 
sites, and therefore J decimations should also be per- 
formed. 

In the following, we concentrate on the active sites, 
and define to each in the log-energy space a range, 

r t = - ln/ij. (6) 

In the improved algorithm we compare the ranges of the 
active sites with the true distances measured between 
them. In this respect two possibilities may happen. 

• If the true distance between two active sites {i and 
j) is not greater than any of their ranges, i.e. <5y < 
r, and (5y < ?y , then i and j are fused together into 
an effective active site, which has a range 

f = n+rj-Sij. (7) 

Also the distance measured from this effective site 
to another site, say k, is given by min(<iifc, djk)- 

• If the range, rj, of an active site, i is shorter 
than any of its true distances from active sites, 
Ti < 5iy, Vj, then this site can not be fused to- 
gether with any other active sites, therefore it is 
turned to 'inactive'. Then we set its weight Z, = 1 
and update the distances, dy,Vj, in Eq.([S]). 

In the above renormalization steps, which can be used 
in arbitrary order, the number of active sites in the sys- 
tem is reduced by one. Repeating these decimation rules 
we arrive to a system having only inactive sites and no 
further fusion steps take place. The complete cluster- 
structure including the excitation energies are readily 
encoded in this configuration, which can be extracted 
without further renormalization steps. In any case the 
final result of the improved algorithm is identical to that 
obtained by the naive SDRG algorithm. 



FIG. 1: (Color online) Snapshots of the naive (left) and im- 
proved (right) SDRG algorithms for the same critical 20 x 20 
sample with box-h disorder and open boundary conditions, 
where 40% of the spins is decimated out. While in the naive 
algorithm numerous new couplings are generated in the im- 
proved algorithm we merely delete sites. 
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FIG. 2: (Color online) Computational time of the algorithm, 
t, as a function of the size of the hypercubic clusters, N, in 
a log- log scale for 2D, 3D and 4D and for the two different 
disorders (fixed-/i +, box-h □). The theoretical prediction, 
t ~ TVlogiV, is indicated by a dashed line. 



In practice one should measure the true distances be- 
tween the active sites. This can be done from a selected 
reference site from which we measure the distances by 
Dijkstra's method^ including one nearby site after the 
other until the range, rj , is reached. While repeating the 
exploration of the shortest paths from all active sites, a 
given site could be crossed from several directions. How- 
ever, it can be shown, that it is always sufficient to cross 
it from one direction onlj«2£. As a consequence, we can 
delete those sites, which have already been explored, from 
the system, because these are no longer needed for the 
calculation. Therefore in a more efficient algorithm ini- 
tially we consider all active sites and perform the mea- 
surements simultaneously and successively delete the ex- 
plored sites. The implementation of this algorithm has a 
time complexity of 0(N log N + E) on any graphs with 
N sites and E edges. A further speed gain is achieved 
by recognizing, that in this parallel implementation the 
paths only need to be explored until reaching a length of 
rj/2 instead of the full range. In Figfljwe compare the 
topology of the renormalized systems in the naive and 
improved SDRG algorithms. 



III. CRITICAL BEHAVIOUR IN DIFFERENT 
DIMENSIONS 

A. Analytical results in ID - a reminder 

In ID the position of the critical point is given by the 
condition^: 5 = [ln/i] av — [lnJ] av = 0, where [. ..] av 
stands for averaging over quenched disorder. At the crit- 
ical point the energy-scale, Q, and the length-scale, L, 
which is the linear size of the system are related as: 
log(flo/f2) ~ > with tp = 1/2. This type of unusual 
dynamical scaling relation is a clear signal of infinite 
disorder scaling. In the paramagnetic phase, S > 0, 
the spin clusters have a finite extent, £, which in the 
vicinity of the critical point diverges as £ ~ with 
v = 2. At the critical point the largest cluster is a frac- 
tal, having a moment: /i ~ L df , with a fractal dimen- 
sion: df = (1 + \/5)/4. From this the magnetization 
exponent is expressed by j3 = xv, with x = d — df and 
d = 1. These SDRG results^, which have been extended 
to the dynamical properties in the Griffiths phases^ have 
been tested by independent analytica l 31 ' 32 and numerical 
calculations 3 ^— . 



We have checked the computational time of the algo- 
rithm, t, for 2D, 3D and 4D hypercubic clusters consist- 
ing of N sites. The results are shown in Fig J2] for the two 
types of randomness. For a given N and for a given type 
of disorder the computation time is practically indepen- 
dent of the topology of the cluster and t is well described 
by the theoretical bound: ~ A log A. (In the studied 
cases the number of edges are proportional to A.) Gen- 
erally, for a given A the renormalization for the box-h 
randomness is faster. 



B. Numerical results for D = 3 and 4 

We have studied finite systems of hypercubic lattices 
in dimensions D = 3 and 4, the largest linear sizes of the 
samples were L — 128 and 48, respectively. (In 2D the 
previously performed numerical investigations^ went up 
to L = 2048.) For each sizes we have renormalized typi- 
cally 40000 random samples (for each type of disorder), 
but even for the largest systems we have treated at least 
10000 realizations. 
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Finite-size critical points 



The precise identification of the critical point of a dis- 
ordered system is a very important issue, since the ac- 
curacy of the determination of the critical exponents de- 
pends very much on it. This question is deeply related to 
the problem of finite-size scaling in random systems 3 ^—, 
since estimates for the critical points are generally calcu- 
lated in finite samples. As has been observed recentl y 37 ' 38 
the key issue in this point is the scaling behavior of the 
finite-size critical points, which are calculated for a large 
set of samples. For a given sample, denoted by a, this 
pseudo-critical point, 8 c (L 7 a), is located at the point, 
where some physical parameter of the system has a max- 
imal value. For example in a classical random magnet 
the susceptibility can be used for this purpose, which is 
divergent at the true critical point of the infinite system. 

For a random quantum system the susceptibility is di- 
vergent in a whole region, in the so called Griffiths phase, 
thus can not be used to monitor finite-size critical points. 
Then in ID random quantum systems the average entan- 
glement entropy turned out to be a convenient quantity, 
the maximum of which is a good indicator of finite-size 
criticality^. In more complicated topology, for ladders^ 
and in 2D systems^ 7 , the so called doubling method^ is 
found to provide an appropriate definition of 8 c (L,a). 
In this procedure one considers two identical copies of 
a given sample, a, which are joined together by surface 
couplings and this replicated sample is denoted by 2a. 
Using the SDRG method one calculates some physical 
quantity (magnetization or gap) in the original and in the 
replicated sample, which are denoted by /(a) and /(2a), 
respectively, and study their ratio, r{a) = f(2a)/f(a), 
as a function of the control parameter, 8. At 8 = 8 c (a, L) 
this ratio has a sudden jump, which is identified with the 
pseudo-critical point of the sample. It has been realized^ 7 - 
that this jump in the ratio is related to a sudden change in 
the cluster structure which is generated during the SDRG 
procedure. For weak quantum fluctuations, 8 < 8 c (L,a), 
between the replicas correlations are generated during 
rcnormalization, which are manifested by the presence of 
a so called correlation cluster. This contains equivalent 
sites in the two replicas. On the contrary for stronger 
quantum fluctuations, 8 > 8 c (L,a), the two replicas are 
renormalized independently For 8 < 8 c (L,a) the mass 
of the correlation cluster, fi(L,a,9), is a monotonously 
decreasing function of 8. Then we identify 8 C (L, a) as the 
point where the correlation cluster disappears. 

Using this doubling method we have calculated pseudo- 
critical points in different dimensions and for different 
forms of the disorder. We illustrate the distributions of 
the pseudo-critical points in FigO for D = 4. To ana- 
lyze these distributions we make use results of finite-size 
scaling theor y 37 ' 38 , which makes statements about the 
average value, 8 C (L), and the width of the distribution , 
A8 C (L), in finite systems of linear size, L. The average 
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FIG. 3: (Color online) Distribution of the pseudo-critical 
points, 9 C (L), for various sizes for fixed-/i randomness for 4D. 
In the inset the scaled distributions are shown as a function 
of y = (6 C {L) - e c )L 1/u ', see the text. 



value of the distribution is expected to scale as: 



(8) 



where 8 C is the true critical point of the system and v s 
is the so called shift exponent. On the other hand the 
width of the distribution is expected to scale as: 



A8 C (L) ~ L- 1 '^ 



(9) 



with a width exponent, v w . We have checked that these 
finite-size scaling relations are satisfied in all dimensions, 
and these relations are used - by comparing results at 
two different sizes (at L and L/2) - to obtain finite-size 
estimates for the exponents. These are presented in Fig|3] 
for the two different forms of disorder. 
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FIG. 4: (Color online) Finite-size estimates for the shift, v a , 
(left) and the width, v w , (right) critical exponents for 2D, 3D 
and 4D and for the two different disorders (fixed-/i +, box-h 
□). The estimated values, as given in Table (I] are indicated 
at the right edge of the figures. 

For a given dimension the estimated exponents are 
found to be independent of the form of the disorder, fur- 
thermore the shift and the width exponents are identi- 
cal within the error of the calculation. Our estimates 
about the critical exponents arc collected in Table]]] to- 
gether with the estimates of the true critical points. 
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TABLE I: Critical properties of the RTIM in different di- 
mensions. In fD the analytical results are from±2 in 2D the 
numerical results are taken fromi£. 





ID 


2D 


3D 


4D 







1.6784(1) 


2.5305(10) 


3.110(5) 




-1. 


-0.17034(2) 


-0.07627(2) 


-0.04698(10) 




2. 


1.24(2) 


0.97(5) 


0.825(40) 


Vs 




1.25(3) 


0.987(17) 


0.74(4) 


X 


3-V5 
4 


0.982(15) 


1.840(15) 


2.72(12) 


V> 


1/2 


0.48(2) 


0.46(2) 


0.46(2) 



We note, that the relation v a = v w = v is character- 
istic for scaling at a conventional random fixed point, 
and the distributions of the pseudo-critical points can 
be rescaled to a master curve in terms of the variable, 
V = {®c{L) ~ O^L 1 /", which is shown in the inset of 
FigH 



dimension which is related to the anomalous dimension 
of the magnetization as i = j3/v = d — df. Compar- 
ing the average mass of the largest clusters at two finite 
sizes, we have calculated effective, size-dependent fractal 
dimensions, as well as effective magnetization scaling di- 
mensions. These are presented in the inset of Fig|5] for 
the different dimensions and using the two disorder dis- 
tributions in Eqs.© and (@|. Extrapolating these values 
the obtained exponents, for a given dimension, do not 
depend on the form of the disorder. These are presented 
in Tableffl 

We have also studied the distribution function of the 
mass of the clusters, Pj,(fi), which is expected to behave 
at the critical point as: Pl(v) = L f P(fj,L f). Accord- 
ing to scaling theory^ P(u) for large arguments has a 

power-law tail, P(u) ~ u~ T , with t = H . In Fig(5]we 

d f 

have plotted P L (fi) for D = 2, 3 and 4 for L = 1024, 128 
and 48, respectively, and good agreement with scaling 
theory is found. 




FIG. 5: (Color online) Left panel: correlation clusters at 
the critical point for fixed-/i randomness in 2D (L = 32) and 
3D (L — 16). Right panel: the connected subgraphs, which 
contain the correlation clusters (see text). 



2. Magnetization 

The magnetization of the system, m(L, 9), is related to 
the average mass of the largest effective clusters, ~p(L, 9), 
as m(L,6) = Jl(L,9)/L d . In the thermodynamic limit, 
L —> oo, in the ferromagnetic phase, 8 = 9 — 9 C < the 
magnetization is finite and vanishes at the critical point 
as: liniL-^oo m(L, 9) ^ [~o~)^ where (3 is the magnetiza- 
tion exponent. At the critical point the largest connected 
clusters arc fractals, which arc illustrated in the left panel 
of FigJS] for 2D and 3D. The mass of these critical clus- 
ters scales as ji ~ L df , where df is the appropriate fractal 
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FIG. 6: (Color online) Distribution of the mass of the clus- 
ters, P L (n), at the critical point for 2D (L = 1024), 3D 
(L — 128) and 4D (L = 48) with box-/i randomness in a 
log-log scale. The scaling results about the asymptotic slopes 
of the curves are indicated by straight lines. Inset: Finite-size 
estimates for the fractal dimension of the critical correlation 
cluster for 2D, 3D and 4D and for the two different disor- 
ders (fixed-/i +, box-/i □). The estimated values, as given in 
Table(T] are indicated at the right edge of the figure. 

Close to the critical point the finite-size magnetizations 
are shown in FigJT] for 3D and 4D. For large L in the 
ordered phases (<5 < 0) the magnetization approaches 
a finite limiting value, whereas for S > it tends to 
zero. In the vicinity of the critical point the finite-size 
magnatizations can be transformed to a master curve, if 
one considers the scaled magnetization, rh = mL x , as 
a function of the scaling variable, 8 = 8L X I V . This is 
illustrated in the insets of FigjT] where the exponents x 
and v are taken from TableJH 
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FIG. 7: (Color online) Finite-size magnetizations in the vicin- 
ity of the critical point for the 3D (left) and the 4D (right) 
models for box-h randomness. In the insets the scaled mag- 
netizations are presented, see the text. 



3. Dynamical scaling 

The dynamical behavior of the RTIM is related to 
the low-energy excitations of the system, the energy of 
which in the SDRG method is given by the values of 
the effective transverse fields at which a given cluster is 
eliminated. Generally, for each such eliminated cluster 
one can define a connected subgraph, which contains the 
given cluster and renormalization of this subgraph gives 
the same energy value. The form of the connected sub- 
graphs of the correlation clusters is illustrated 4 ^ in the 
right panel of Fig|5] The energy-parameter of a given 
sample, which is denoted by e(L, a), is given by the small- 
est effective transverse field, not considering the trans- 
verse field of the correlation cluster, if it exists in the 
system. The distribution of the log-energy parameter, 
j(L, a) = — loge(L, a), is shown in the upper panel of 
Fig |S] at the critical points of the 3D and 4D systems. The 
distributions for both dimensions are broadening with in- 
creasing L, which is a clear signature of infinite disorder 
scaling. The typical value of the log-energy parameter 
grows with the size as 'y(L) ~ L^' , thus the appropriate 
scaling combination is given by: 7 = (j(L) — jo)L~^. 
Here ip is a scaling exponent and 70 is a non-universal 
constant. The scaled distributions are shown in the lower 
panel of Fig|5] 

We have calculated effective, size-dependent ip expo- 
nents, by comparing the widths of the distributions of the 
log-energy parameters at two sizes (L and L/2). These 
are given in Fig |H] for the different dimensions and for the 
two different form of disorder. As for other exponents 
studied before, the estimates for tp for a given dimension 
do not depend on the actual form of the randomness. 
These are summarized in Tablcjll Interestingly the ex- 
ponents for all studied dimensions are close to 1/2, which 
is the exact value in ID. This observation can be ex- 
plained by the fact, that the connected subgraphs, which 
are related to the energy-parameter of the sample, are ba- 
sically one-dimensional objects in all studied dimensions, 
see in the right panel of FigJS] The very small variation 
of ip with the dimensionality is probably due to the fact 
that the renormalizcd couplings and transverse fields of 
the connected subgraphs are more and more correlated 
in higher dimensions. 
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FIG. 8: (Color online) Distribution of the log-energy param- 
eters as the system size in 3D and 4D for box-h randomness 
at the critical point (upper panel). The scaled distributions 
are shown in the lower panel, using the non-universal param- 
eters: 70 = 3.1 (3D) and 70 = 4.06 (4D), as described in the 
text. 



Finally we note that thermodynamic singularities at 
a small temperature, T, but 5 = are related to the 
critical exponents in Table|I] For example the suscep- 
tibility, x, and the specific heat, Cy, behave a o 11 ' 14 : 
X(T) ~ (logT)( d - 2a; )/^/r and C V (T) ~ (logT)"^. 
The similar relations at T = but with a small longi- 
tudinal field, H, arc given by: X (H) ~ (\ogH)- x ^/H 
and C V (H) ~ (log ifT^. 




FIG. 9: (Color online) Finite-size estimates for the if) critical 
exponent for 2D, 3D and 4D obtained from the width of the 
log-energy distributions for the two different disorders (fixed- 
h +, box-h □). 
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4- Griffiths effects 

We have also studied the distribution of the low-energy 
excitations outside the critical point. In the paramag- 
netic phase the distribution of j(L) for different sizes 
are shown in FigQJJ and in FigfTTl for the 3D model 
(S = 0.37) and for the 4D model (6 = 0.34), respectively. 
As seen in these figures the distributions have approxi- 
mately the same width, they are merely shifted with in- 
creasing L. This behaviour is in agreement with scaling 
theory in the disordered Griffiths-phase^, where the typ- 
ical excitation energy scales with the size as: e(L) ~ L~ z , 
where z = z(S) is the dynamical exponent, which depends 
on the distance from the critical point, S. Then the ap- 
propriate scaling combination is: 7 = j(L) — z log L — 70, 
in terms of which a scaling collapse of the distributions 
are found, which is shown in the inset of FigfTU] and 
FigfTT] One way to estimate z(6) is to analyze the scaling 
collapse of the distributions, or equivalently to compare 
the shift of the distributions with L. From this type of 
analysis we obtain d/z = 0.75(4) for the 3D model and 
d/z = 0.80(4) for 4D. There is, however, another possi- 
bility to calculate d/z from the asymptotic form of the 
distributions. If the low-energy excitations are localised, 
which is satisfied for the RTIM, the distribution function 
of the scaled variable, 7, is expected to follow extreme 
value statistics^ and given by the Frechet distribution 4 ^: 

lnp(7-7o) = -~7-exp + in ■ ( 10 ) 

Indeed the scaled distributions in the insets of FigslTOl 
and [TT] are well described by this form, having just one 
free parameter, j . From Eq. ([T0"]) follows that the asymp- 
totic slope of lnp(7) vs. 7 is just d/z, what we have 
measured in FigsfTU] and [TT] The estimates for d/z for 
different values of L, as given in the captions are in good 
agreement with our previous estimates from the shift of 
the curves. 

We have repeated these type of calculations at other 
points of the disordered Griffiths phase and we have 
measured i5 dependent dynamical exponents. We could 
not, however, check the scaling resul t 11 ' 14 for small 5: 
d/z cx 5 U ^' , due to strong finite-size effects in the vicinity 
of the critical point. (In 2D this type of analysis has been 
performed ir>i£.) 

We have also calculated the distribution of the log- 
energy parameters in the ordered Griffiths phase, which 
is illustrated in FigfT2]for the 3D model at 5 = —0.33. In 
the ordered phase there is a huge magnetization cluster, 
which has a very small effective field and the energy pa- 
rameter is given by the second smallest effective field of 
the RG process. As in the disordered Griffiths phase the 
width of the distributions is approximately L indepen- 
dent and the distributions are shifted with L. However, 
the amount of shift with L, as well as the shape of the dis- 
tributions arc different in the two cases. This is related 
to the scaling result, that the typical value of the exci- 
tation energy in the ordered Griffiths phase scales with 




FIG. 10: (Color online) Distribution of the log-excitation 
energies in the disordered Griffiths-phase of the 3D model 
for box-/i disorder at 8 = 0.37 in a log-lin scale for different 
sizes. The slopes of the straight lines indicating the tail of the 
curves are d/z = 0.71(2), 0.69(2) and 0.71(2), for L = 24,48 
and 96, respectively. In the inset the scaled distributions are 
shown with d/z = 0.73, which is well described by the Frechet 
distribution with 70 = 6.09 (full line). 




FIG. 11: (Color online) The same as in FigfJO] for 4D at 
S = 0.34. The slopes of the straight lines indicating the tail of 
the curves are d/z = 0.81(3), 0.81(3) and 0.76(3), for L = 8, 16 
and 32, respectively. In the inset the scaled distributions are 
shown with d/z — 0.8, and the full line represents the Frechet 
distribution with 70 = 5.60. 



the size as: e(L) ~ — In ' (L), thus the appropriate scal- 
ing combination is: 7 = j(L) — Aln 1 ^(L) — 70, with A 
and 70 being nonuniversal constants. Using this variable 
the distributions have a scaling collapse as shown in the 
inset of FigIT2l We should, however, mention that due 
to finite-size effects we can not obtain an independent 
estimate of the exponent of the logarithm, being theoret- 
ically 1/d. We are facing to the same kind of limitations 
concerning the shape of the scaling curve, which asymp- 
totically should behave as: \np(/y) ~ — 7 d , according to 
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scaling theory. Our data, however, are still not in the 
asymptotic regime. 
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FIG. 12: (Color online) The same as in Fig[l0]for 3D at the 
ordered Griffiths phase, S — —0.33. In the inset the scaled 
distributions are shown with A = 10.0. 

Finally we note that in the disordered Griffiths phase 
the singularities of the susceptibility and the specific 
heat at a small temperature are given b y 11 ' 16 : x(^") ~ 
T -i+d/z and c v (T) ~ T d l z . The same expressions in the 
ordered Griffiths phase are: x(T) ~ exp(-C| logT| d )/T 
and Cv(T) ~ exp(-C"| logT| d ). 



C. Numerical results for Erdos-Renyi random 
graphs 

Here we come back to the question posed in the In- 
troduction about the possible value of the upper critical 
dimension, D u , in the problem. In order to answer to this 
question we consider Erdos-Renyi (ER) random graphs^ 
with a finite coordination number, which are represent- 
ing the large-dimensional limit of our lattices. Generally 
an ER random graph consists of N sites and kN/2 edges, 
which are put in random positions. In order to have a 
percolating random graph we should have k > 1. Here 
we have used k = 3, but some controlling calculations 
had also been done with k = 4. In the actual calcula- 
tion we have put the RTIM on ER random graphs and 
study their critical behaviour by our improved algorithm 
of the SDRG method. Due to infinite dimensionality of 
ER clusters we had to modify some parts of the analysis 
used in Scc lIII Bl for finite D. 

As for D < 4 we have calculated sample dependent 
pseudo-critical points, but now in the doubling method 
the two identical copies of the sample have been con- 
nected by N/2 random links. The distribution of the 
calculated pseudo-critical points for different values of N 
are shown in Fig|13j The general behaviour of the distri- 
butions is similar to that for finite-D, see Fig J3] for 4D, 
but in the present case for large 6 C values there is an N- 
independent background of the distributions. This back- 
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FIG. 13: (Color online) Distribution of the pseudo-critical 
points, 6 C (N), for the ER random graphs with box-h ran- 
domness. In the inset the scaled distributions are shown as a 
function of y = (6 C {N) - 6 c )N 1/u with u = 6. 



ground is probably due to the large number of connect- 
ing random links between the replicas. This background, 
however, has a very small weight to the distributions and 
does not influence the analysis of the properties of the 
pseudo-critical points. Concerning Fig ll3l we have mea- 
sured the shift, \0 c -T c (N)\ ~ jV- 1 /^ an d the width 
A8 C (N) ~ JV _1 /" lu of the distributions, in analogy with 
the finite dimensional problem. (Compare with Eqs.® 
and as well as with v s — > uj s /d and v w u> w /d, re- 
spectively.) From two-point fits we have calculated effec- 
tive exponents (see Fig ll4[) from which we have obtained 
the estimates, uj s = 4.5(1.5) and uj w = 7.8(2.0), which are 
valid for both type of randomnesses. We note that the 
relative error of the estimates is somewhat larger than for 
finite D, but still the two exponents of the distribution 
agree with each other giving u> = 6.(2). Using the scaled 
variable, y = (6 C (N) — B^N 1 ^ , the distributions show a 
scaling collapse as illustrated in the inset of FiglT^l 

We have studied the fractal properties of the correla- 
tion cluster, the average mass of which is found to scale 
at the critical point as: ~p{N) ~ iV 5 . From two-point fit 
we have obtained effective values for which are shown 
in FigfMland which are extrapolated to $ = 0.17(5). We 
note that according to scaling theory the magnetization 
exponent of the RTIM on the ER random graph is given 
by: /3 = u;(l-i?) =5.(2). 

We have also investigated the distribution of the log- 
energy parameters at the critical point, which is shown 
in Fig[l5] for different sizes of the ER clusters. In or- 
der to see the possible existence of infinite disorder scal- 
ing we have measured the width of the distributions, 
which are shown in the inset of Fig |15l as a function 
of logiV. As seen in the inset the width of the distri- 
bution can be parametrized as Wo + W\ log 6 N, where 
the constant is Wo ss for fix-/i randomness and it is 
Wq ~ 1.2 for box-ft, randomness. In both cases the ex- 
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FIG. 14: (Color online) Finite-size estimates for the critical 
exponents in ER random graphs for the two different disorders 
(fixed-/i +, box-h □). The estimated values for large- N are 
indicated at the right part of the figure. 
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FIG. 15: (Color online) Distribution of the log-energy pa- 
rameters as a function of the size of the ER clusters for box-h 
randomness at the critical point. In the inset the width of 
the distribution is shown as a function of log TV for the two 
different disorders (fixed-/i +, box-h □). 



poncnt in the logarithm can be estimated as: e = 1.3(2), 
thus the increase of the width of the distribution is some- 
what larger than linear in log TV. This fact justifies that 
even for ER random graphs the critical behaviour of the 
RTIM is controlled by a (logarithmically) infinite disor- 
der fixed point^. Thus we conclude, that the upper crit- 
ical dimension of infinite disorder scaling of the RTIM 
is D u = oo. For ER random graphs the singularities of 
the susceptibility and the specific heat at small temper- 
ature are given by: Iog[Tx(T)] - (2d - 1)| logTl 1 / 6 and 
loglCviT^-llogTl 1 ^. 



IV. DISCUSSION 

In this paper we have considered the random 
transverse-field Ising model, which is a basic model of 
random quantum magnets and studied its critical be- 
havior in different dimensions by a variant of the SDRG 
method. These investigations are made possible that we 
have developed an improved numerical algorithm for the 
SDRG method so that we could renormalize clusters with 
up to TV ~ 4 x 10 6 sites irrespective of their dimension- 
ality and topology. We have found strong numerical ev- 
idence that the critical behaviour of the RTIM for all 
dimensions up to D = 4 is governed by infinite disorder 
fixed points. This fact justifies the validity of the use of 
the SDRG method as well as indicates that the obtained 
critical properties of the model, which are summarized in 
Table|T]are asymptotically exact. This means that with 
increasing sizes in the calculation the critical exponents 
approach their exact value. We have demonstrated by 
using different disorder distributions in the initial mod- 
els that the strong disorder fixed points are universal, the 
critical parameters do not depend on the actual form of 
the disorder. We have also studied the behaviour of the 
systems in the vicinity of the critical points and good 
agreement with scaling considerations are obtained. 

We have considered the upper critical dimension of infi- 
nite disorder scaling of the RTIM and studied the critical 
behaviour of the model on Erdos-Renyi random graphs 
by the improved SDRG algorithm. Our results indicate 
that even in this, formally infinite dimensional lattice the 
critical behaviour of the RTIM is governed by a (loga- 
rithmically) infinite disorder fixed point, thus the upper 
critical dimension is D u = oo. 

Our results presented in this paper are relevant for sev- 
eral other problems, too, since the IDFP of the RTIM is 
expected to govern the critical properties of a large class 
of random systems, at least for strong enough disorder. 
These are, among others, random quantum ferromagnetic 
systems having a continuous phase transition at which a 
discrete symmetry of a non-conserved order parameter 
is broken. Examples are the quantum Potts and clock 
models^ as well as the Ashkin- Teller model^. Also the 
quantum spin glass (QSG) problem could be related to 
the IDFP of the RTIM. For the QSG the distribution 
of couplings in Eq.@ contains antifcrromagnetic terms, 
too, however, at an IDFP frustration is expected to be 
irrelevant. Thus the critical exponents in Tabled with 
some appropriate modifications of the scaling relations in 
the ordered phased should hold for the QSG, at least for 
strong enough disorder. Also noncquilibrium phase tran- 
sitions in the presence of quenched disorder are expected 
to belong to the universality class of the RTIM^i and the 
random walk in a self-affine random potential^ could be 
related to the RTIM. 

Finally we mention that the ideas about the numeri- 
cal implementation of the SDRG method in Scc lII Al can 
be generalized for another models. In the Appendix we 
outline the elements of the improved SDRG algorithm 
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for the random quantum Potts model^, as well as for 
the model of disordered Josephson junctions^. These 
results can be used to investigate the critical behaviour 
of these systems in higher dimensions. 

Appendix A: Improved SDRG algorithm for other 
models 

The SDRG approach has been applied for a series of 
random quantum and classical problems mainly in one 
dimension. In higher dimensions the numerical imple- 
mentation of the SDRG method for these models has ba- 
sically the same problems as the naive algorithm for the 
RTIM. In these cases one can try to apply and generalize 
the concept of our improved algorithm. Here we present 
the appropriate RG rules for two random quantum mod- 
els: for the disordered g-state quantum Potts model and 
for the disordered quantum rotor model, which is a stan- 
dard model of granular superconductors and Josephson 
arrays. 

1. Disordered q-state quantum Potts model 

In this model at each lattice site, i (or j) there is a q- 
state spin variable: Si = 1,2, ... ,q and the Hamiltonian 
is given byi&: 

h = - E J ^ s ^ E - E M " ■ ( A1 ) 

(ij) i q k=l 

Here the first term represents the interaction between the 
spins and the second term is a generalized transverse field 
where Mi is a spin- flip operator at site i: M,|sj) = |si + l, 
mod q). As for the random transverse-field Ising model, 
what we recover for q = 2, the Jy couplings and the 
hi transverse fields are random variables. The SDRG 
decimation rules are very similar to that of the RTIM as 
described in Sec|Hl which differ only in an extra factor: 
k = 2/q. 

J-dccimation : the effective transverse fields are given 
by: h = nh i hj/J i j. 

fa-decimation : the effective couplings arc given by: 
Jjk = KJjiJik/hi, which should be supplemented by the 
maximum rule. 

In the improved algorithm in Scc lII Al the distances 
and ranges in the log-energy space, see Eqs.dS]),© and 
(J7J), are extended by a constant: do = ln(g/2) = — ln/c, 
which now read as: 

dij = - In J i:j + | In hi + | In hj + d , ( A2) 

n = — ln/ij + d , (A3) 



f = r i + r j -Sij+do. (A4) 

This generalization works for 2 < q < oo. 

2. Disordered Josephson junctions 

Here we consider disordered bosons with an occupation 
operator, hi, and a phase- variable, (pi, at site i. The 
system is described by the following Hamiltonian^: 

H = ~E Jij cos (ifi - <pj) + ^2Uih? , (A5) 

(ij) « 

with random Josephson couplings and Ui charging 
energies. The SDRG approach has also been applied to 
this model resulting in the following RG rules. 

t/-dccimation: If the strongest parameter in the Hamil- 
tonian is a grain charging energy Ui, then site i is elimi- 
nated and effective couplings are generated between the 
nearest neighbours of i, say j and k. In second-order per- 
turbation calculation this is given by: Jjk = JjiJik/Ui, 
which have to be supplemented with the maximum rule. 

J-decimation: If the strongest coupling in the system 
is a Josephson coupling, Jy, then the two sites form a 
composite site having an effective charging energy, Ui, 
which does not depend on the value of Jy but given by: 
111 

Wi ~ ul + U]' 

In the improved SDRG algorithm the distance and the 
range in Eqs.([5]) and ([6]) are modified with the substitu- 
tion hi — > Ui as: 

dij = - In Jij + | In Ui + | In U 3 , ( A6) 

n = -]nU i . (A7) 

On the contrary the updated range in Eq.([7|) has a dif- 
ferent form and given by: 

f = In [cxp (rj) + exp (fj)]. (A8) 
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